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I. INTRODUCTION 

l> : 

Q\ , The one-loop corrections of quantum electrodynamics introduce nonlinearities in the equations of the electromag- 
0^ ' netic field. These corrections manifest themselves through an index of refraction, electric permittivity and magnetic 
permeability tensors which are a function of field strength Jl],|| . The vacuum responds to an applied field like a non- 
linear paramagnetic substance || . Using an analytic expression for the effective Lagrangian of QED to one- loop order 
Q , we derive the magnetic magnetic permeability tensor as a function of the applied magnetic field and calculate the 
^! . one-loop corrections to the field of a macroscopic magnetic dipole. 

o : 

. II. THE EFFECTS OF NON-LINEARITY ON THE FIELD EQUATIONS 

' The effective Lagrangian of QED at one-loop order consists of the sum of a linear and a non-linear term 

00 : 

VO ■ C = C +Cx- (1) 

m 

Both terms of the Lagrangian can be written in terms of the Lorentz invariants, 

o 

l> ' 1= F^F^ = 2 (|B| 2 - |E| 2 ) (2) 

G\ ; 

"*^h , and 

£ K = -(±e^F Xp F,}j = -(4E-B) 2 . (3) 

'-^ ' following Heisenberg and Euler ||. The Greek indices count over space and time components (0, 1, 2, 3). 

. Since we are interested in the static properties of a magnetic field, we can take K = 0. If we apply the Euler-Lagrange 
' condition to extremizc the action, we obtain, 

: V x H = (4) 

H = -4^B (5) 

where the factor of —4 is inserted for later convenience. 

Since V x H = 0, we take H = — V</>, where <f> is the magnetic scalar potential. Furthermore, the field B is derived 
from a vector potential (i.e. B = V x A) and so we also have, 

V ■ B = 0. (6) 

If the relationship between H and B were linear, this equation would be satisfied by V 2 </> = 0. However, we will 
assume a small non-linearity between the two fields, 

-4^-^ + M!(B 2 ) (7) 

where fi' is a constant and fi[ a function such that fi[(B 2 ) fi' . 
We can invert the relationship between the two fields to first order 

B= ( Al0 + m (H 2 ))H (8) 
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Now we recast the field equation with the magnetic potential 

V- [( Mo + m (tf 2 ))V0] =0, 
write cj) = (f>o + 4>i j and solve the equation order by order 
V 2 </) = 

V 2 0i = peff = -V 

where 
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For a magnetic dipole, 



0o (r) 



m • r 



V^o • (V0 O • V)V0 o = 3 



5(m-r) z + 3|m| 2 |r |2 



m • r 
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(9) 

(10) 
(11) 

(12) 

(13) 
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or more conveniently in spherical coordinates where we have taken the dipole moment m to be aligned along the 
z— axis, 



V S r z 



V0 O • (V0 O • V)V0 o = 12^^ 



101 



(15) 
(16) 



III. THE LAGRANGIAN TO ONE-LOOP ORDER 

Heisenberg and Euler j|] and Weisskopf [pjindependently derived the effective Lagrangian of the electromagnetic 
field using electron-hole theory. Schwinger [£[ later rederived the same result using quantum electrodynamics. In 
Heaviside-Lorentz units, the Lagrangian is given by 



4) = -^ (17) 




where B k = E k = sa 2.2 x 10 15 Vcm _1 w 4.4 x 10 13 G. Dittrich and Reuter 0] have derived the second-order 
corrections to the Lagrangian and found them to be in general an order of a smaller than the one-loop corrections 
regardless of field strength; consequently, the one-loop correction should be adequate for all but the most precise 
analyses. 

In the weak field limit Heisenberg and Euler || give 



We define a dimensionless parameter £ to characterize the field strength 

and use the analytic expression of this Lagrangian for K = derived by Heyl and Hernquist 
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where 



X {x) = 4 
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ln(I> + l))dv + - In 



21n47r- (4 In A + - In 2) 
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(22) 



where 



lnA = — - C (1) (-l)- 
12 ^ V ' 



(23) 



Here (a;) denotes the first derivative of the Riemann Zeta function. 

With the analytic form of the Lagrangian, calculating fi' and p! x is straightforward and a provides a convenient 
ordering parameter. 
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where £ = B/Bk, and 



(1) _ dXg(x) 

[X) ~ dx ■ 



Inverting this relationship to first order yields, 
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£ x (i) ( Bk 
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where ip(x) is the digamma function, 



tp(x) 



dlnT(x) 
dx 
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(29) 

(30) 

(31) 



The expression for /zi agrees numerically with the results of Mielnieczuk et al. ||. The function fi\ (if 2 ) may 
conveniently be expanded in the weak-field limit (H < Bk/2), 



£\H 2 



a 8 ^ 2 2 ^ 2(j+2 ) f H ^ 2j 
2nB 2 ^ 2j + 3 



5* 



(32) 



where Bj denotes the jth Bernoulli number, and in the strong-field limit (if > Bk/2) 
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Q 1 J lB 'l lB l 



(H 1 = — — — / « _ fc 
1 ' 2tt 1 3 H 2 2H 3 



1 — ln7r 
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i=5 



2J-3 7 



:C(J 



(33) 



where we have used the expansions of Rcf . [Q . 

As apparent from Fig. [TJ, p^'(H 2 ) is constant up to approximately H = 0.5Bk and then begins to decrease quickly 
as H~ 2 . The existence of these two regimes allows us to find analytic solutions for the correction to the potential 



IV. SOLVING FOR THE FIRST-ORDER CORRECTION 
A. Weak-field limit 

p^(H 2 ) is constant as long as H <C B^. In this regime Eq. [TT|may be solved analytically. Since spherical harmonics 
are eigenfunctions of the angular component of the Laplacian operator, it is expedient to expand the right-hand side 
of Eq. [y] in terms of spherical harmonics |j) . 

Plm(r)= smddO / dcj ) Y l * m (e, ( j ) )p cS (r,6, ( f>) (34) 
Jo Jo 

With this definition the correction to the potential is given by M 



l.m 



(35) 



If Pi (H 2 ) is a constant, we see from Eq. |Tl| and Eq. |l6| that the expansion given by Eq. |34| is straightforward 

(1) 3 

p 10 (r)=-24>/3^-^ (36) 
Po r w 

/— (i) s 
7r a \ m 
(37) 

If we combine these results with Eq. we see that the interior integral diverges if p^ is a constant. This does not 
present a problem if we insert an inner bound (ro) to the interior integral. This inner bound is defined as the radius 

at which either p^p (H 2 ) begins to change (i.e. when H > O.bBk) or when the zeroth order potential is no longer 
given by the dipole formula (i.e. at the surface of the object). We obtain 

0i,io(r, 0, 0) = - T 4 " T 0) (38) 



9 /io ?" 2 V^o r6 

7T ^ TO 3 / 1 1 I I 



01,30 (r, M) =6W-^ y 3O (0,0). (39) 

V 7 ^ r 4 V 7 r o 11 r J 

These functions describe radially dependent corrections to the dipole and hexapole moments of the object under 
consideration. If we define the higher moments in analogy to Eq. rL5l, 



we obtain the following corrections 



= \/ 2f + lr ffW g 'fl' (40) 
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To eliminate the dependence on the inner bound, ro, we assume that the dipole and hexapole moments are known at 
a radius r s > r$ and calculate the difference between the known moments at r s and those measured at infinity. 



Substituting the values of /t and (jfy~ , we obtain, 



/ \ / \ a 8 a / m 1 

m(r = 00)- m(r = r s ) = Am = — -m -- 



M 30 (r = 00) - M 30 (r = r.) = AM 30 = A JL mr * g 2- 



(43) 
(44) 



Both the corrections are given in terms of the magnetic field strength at r s . Since we have assumed that fj,^ (H 2 ) 
is constant near r s , the field strength at r s must appreciably be less than B k for this set of approximations to be 
valid. Consequently, the correction to the dipole moment is indeed quite small, less than one-thousandth of the "bare" 
dipole moment. However, there is a correction to the hexapole moment even without a "bare" hexapole moment. 
Thus radiative corrections generate a hexapole field which is in principle measurable at infinity. 



B. The General Case 



In the strong-field limit, we must use the general equation (Eq. |o|) for ^(H 2 ) because the magnetic field strength 
varies as a function of 9 around the dipole. Before tackling this problem numerically, we can glean several character- 
istics of the solutions from Eq. |ll| and Eq. [}4| Since for a dipole p c g is a odd function of 9 and constant with respect 
to <f), the contributions, 



Pim{f) = if / is even or m ^ 0. 
Furthermore, we can more conveniently write 

I (X 1 TV? 
Plm(r) = — ^-^2— Xlm(P) 

fj, 2n Bl r w 

where Xim(P) is a dimensionless function of a dimensionlcss argument, 

' r*B k 

For a dipole, 



(45) 



(46) 



(47) 



X i m (P) = -2 / sm0d0 / d01^(6\. 
Jo Jo 



2tt 



a 1 



__) ( 3 cos 6> (5 cos 2 6» + 3) J/4 ' [P 2 B 2 k (3cos 2 6> + l)] . (48) 



(i) 



We calculate numerically the functions Xlm(0) for the first three odd harmonics and depict the results in Fig. 0. 
The limiting expressions for the weak-field limit are easily calculated, 



Xw(P) = V3tt 
X3o(/3) = V7w 



32 13312^ _ 167936^ + 
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In the strong-field limit, we have the following approximations 



(49) 
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Xw(fi) = - (27V3 - 4^J (3- 2 + O (p~ 



X30 



09) = 243^ ( 81 - U ^ n ) + ° 3 " 8 ) 



xM = - 



1215 



lbr (l863 - 340V37r) /T 2 + O (/3 



(52) 
(53) 
(54) 
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The integrals for each spherical-harmonic component of the first-order correction (Eq. p5| ) may be recast in terms of 
integrals over j3 



h{r,t 



1 TO J__£*o2 Ylm 

"3r 2 /xo27r^ f- 



2/ + 1 



/So 



/? (i - 7)/3 / » m M^ ( '- 4)/3 d W + r 0+8)73 / x im w^ ( ' +5)/3 d W 



(55) 



where the first integral is to be evaluated in the limit as /3q — ► oo. Although Eq. |55| appears to be scale free, the 
cutoff P has a physical interpretation. Firstly, it can be taken to be the magnetic field strength at the surface of the 
object. For a point magnetic dipole (e.g. an electron), the interpretation is more subtle. As one approaches a point 
dipole, not only does the field become arbitrarily strong, so do the field gradients. When these gradients become 
larger in magnitude than B^j\ e (A e = h/m e c, the electron Compton wavelength), the Heisenberg-Euler Lagrangian 
is no longer applicable; therefore, we do not expect our expressions for p c ff to be valid arbitrarily close to a point 
dipole. The radius or field strength at which our expression for p c g fails depends on the intrinsic dipole moment of 
the object (m), 
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As in the weak-field case, we calculate the shift in the observed multipole moments at infinity relative to the known 
moments at some inner radius r s > ro or equivalently at some field strength (3 S < (3q. 
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(60) 



Fig. |^ depicts the shifts in the moments between the surface of the dipole and infinity. 

For the higher moments (I > 1), the integrals in Eq. ^5] are well behaved so we need not set a inner bound (/3q) and 
we can define in general 



M w> i(0) = - 



m \ 



3(2/ + 1) \B k J 



(i-l)/3 



la 21 + 1 



HQ 2TT V 47T J a 
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This integral converges for all / > 1. For Z = 1, it diverges logarithmically as v — > oo. 



V. CONCLUSIONS: APPLICATION TO NEUTRON STARS 



The environment of a neutron star is strongly magnetized; therefore, the one-loop QED corrections may be signif- 
icant, especially for those neutron stars with ultrastrong surface fields (magnetars) |l(]]. Although plasma prob ably 
fills this region, QED vacuum corrections dominate the contribution of the plasma to the magnetic permeability [p_ 1| . 

The structure of the magnetic field at the surface of a neutron star is an important clue to the origin of neutron-star 
magnetic fields. Several authors have proposed [jl2| that currents in the thin crust generate the observed magnetic 
fields. In this case, the field structure will be dominated by high-order multipoles with / ~ 8r c /r* 3> 1 where 8r c is 
the thickness of the crust and r* is the radius of the star . Arons |l3| argues further that the observed spin-up 
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line for millisecond pulsars constrains the strength of higher-order multipoles at the surface to be no more than 40 % 
of surface strength of the dipole. 

The current results complicate this argument. The location of the observed spin- up line and the value of spin-down 
index of a pulsar depend on the strength of the various moments of the magnetic field at the light cylinder. Our 
results show that the vacuum itself may generate higher magnetic moments between the neutron star surface and the 
light cylinder. Fig. [| depicts the fractional contribution of higher magnetic moments to the field strength at the light 
cylinder as a function of the surface dipole field strength and the pulsar period. 

Even for magnetars near their birth, the vacuum adds only a small correction to the field strength at the light 
cone. Because of the weakness of the QED coupling, the one-loop corrections to otherwise classical descriptions of a 
magnetic dipole tend to be small for all but the most extreme field strengths. 

Observing this effect would be difficult. For a magnetar, measurements of the field strength at two different radii 
and an estimate of the strength of higher order moments at the surface each to a precision of one part in one-thousand 
would be required. However, if one could argue that an object had no hexapole or higher-order multipolc intrinsically, 
a measurement of a higher-order multipolc far from the object would uncover the effects of one-loop corrections. For 
example, an electron is intrinsically a magnetic dipole. If one ignores the contribution of terms in the Lagrangian 
depending on field gradients, one would expect the vacuum surrounding an electron to generate higher order multipole 
fields. Determining whether this occurs is beyond the scope of this work. 
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FIG. 1. ^ ] (H 2 ) as a function of £, H 
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FIG. 2. Xlm(P) as a function of (3, H. The solid lines indicate where Xlm(/3) is positive. The dashed lines indicate negative 
values of the ordinate. 
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FIG. 3. The difference in the strength of the multipole moments measured at infinity (/3 — * 0) and at the surface, 
Mio(0) — Mio(/3 3 ) (AM(o), as a function the strength of the dipole field at the surface, /3 3 ,H 3 . The solid lines indicate 
the moments for which the vacuum acts paramagnetically, i.e. AM 10 (/3 S ) > 0. The dashed lines indicate negative values of the 
ordinate. 
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FIG. 4. The fractional contribution of higher magnetic moments to the field strength at the light cylinder as function of 
surface dipole field strength for periods of one millisecond and one second. The contribution of the vacuum decapole at the 
light cone for P = Is is too small to be depicted on this graph. 
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